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ABSTRACT 

We propose a ray-tracing method to estimate gravitational waves (GWs) gen- 
erated by anisotropic neutrino emission in supernova cores. To calculate the 
gravitational waveforms, we derive analytic formulae in a useful form, which are 
applicable also for three-dimensional computations. Pushed by evidence of slow 
rotation prior to core-collapse, we focus on asphericities in neutrino emission and 
matter motions outside the protoneutron star. Based on the two-dimensional 
(2D) models, which mimic SASI-aided neutrino heating explosions, we compute 
the neutrino anisotropics via the ray-tracing method in a post-processing manner 
and calculate the resulting waveforms. For simplicity, neutrino absorption and 
emission by free nucleons, dominant processes outside the PNSs, are only taken 
into account, while the neutrino scattering and the velocity-dependent terms in 
the transport equations are neglected. With these computations, it is found 
that the waveforms exhibit more variety in contrast to the ones previously es- 
timated by the ray- by-ray analysis (e.g., Kotake et al. (2007)). In addition to 
a positively growing feature, which was predicted to determine the total wave 
amplitudes predominantly, the waveforms are shown to exhibit large negative 
growth for some epochs during the growth of SASI. These features are found 
to stem from the excess of neutrino emission in lateral directions, which can be 
precisely captured by the ray-tracing calculation. Reflecting the nature of SASI 
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which grows chaotically with time, there is little systematic dependence of the 
input neutrino luminosities on the maximum wave amplitudes. Due to the nega- 
tive contributions and the neutrino absorptions appropriately taken into account 
by the ray-tracing method, the wave amplitudes become more than one-order-of 
magnitude smaller than the previous estimation, thus making their detections 
very hard for a galactic source. On the other hand, it is pointed out that the 
GW spectrum from matter motions have its peak near ~ 100 Hz, reflecting the 
SASI-induced matter overturns of 0(10) ms. Such a feature could be characteris- 
tic for the SASI-induced supernova explosions. The proposed ray-tracing method 
will be useful for the GW prediction in the first generation of 3D core-collapse 
supernova simulations that do not solve the angle- dependent neutrino transport 
equations as part of the numerical evolution. 

Subject headings: supernovae: collapse — gravitational waves — neutrinos — 
hydrodynamics 



Introduction 



No longer gravitati onal- wave astronomy is a fantasy . In fact, gravitational wave de- 
tectors, such as LIGO (|Thornelll99,4 kbbott et allboO^^h . VIGRoQ, GEO60C@, TAMA300 



( lAndo fc the TAMA collaboration 



20021: 



AndoetaL 



2005h . and AIGol with their interna- 



tional network of the observatories, are beginning to t ake d ata at sensitivities where astro- 
physical events are predicted (see, e.g., iHough et al.l (120051 ) for a recent review). For the 
detectors, core-collapse supernovae especially in our Galaxy, have be en proposed a s one o f 
the most plausible sources of gravitational waves (see, for example, iKotake et al.l (l2006l ): 
OttI (120091 ) for recent reviews). Since the gravitational waves (plus neutrinos) are the only 
tool that gives us the live information of the central engine of core-collapse supernovae, the 
detection is important not only for the direct confirmation of general relativity but also for 
disclosing the supernova physics itself. 

Traditionally, most of the t heoretical predictions of g r avitational waves (GWs) have fo- 



cused o n the bounce signals (e.g..lZwerger fc Miillerl ( 



(120041 ) : l0tt et all (l2004[ ): IObergaulinger et al.l (120061 ): 



997h:lKotake et al. (120031) : Shibata k Sekiguchi 



Cerda-Duran et al.l (120071 ) : lDimmelmeier et al. 



"'^ http://www.ego-gw.it/ 
^http://geo600. aei.mpg.de/ 
ht t p :// WW w . gravity, uwa . edu . au/ 
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( 120071 ): IScheidegger et al.l (120081 ): iDimmelmeier et al.l (120081 ) and references therein). How- 
ever recent stellar evolution calculations suggest that rapid rotation assumed in most of the 



previ ous studies is not canonical for the progenitors with neutron star formation (IHeger et al. 



20051 ). To explain the observed rotation periods of radio pulsars, t he rotation perio ds of the 
iron core before collapse are estimated to be larger than ~ 100 sec (lOtt et al.ll2006l ). In such 
a slowly rotating case, the detection of the bounce signals becomes very hard even by the 
next-generation laser interferometers for a Galactic supernova (e.g.. iKotake et al.l (120041 )). 



Besides the rapid rotation of the core, convective matter motions and anisotropic neu- 
trino emission in the much later postbounce phase are expected to be the primary GW 
sources with comparable amplitudes to the bounce signals. Thus far, various physical in- 
gredients for producing asphericities and the resulting GWs in the postb ounce phase, have 
been studied such as the r o les of pre-c o llapse density inh omogeneities (IBurrows fc Hayes 



1996 



Miiller fc Janka 



19971 : iFryerl |2004J : iFryer et al.ll2004l ). moderate rot ation of the iro n 



core (IMiiller et al.ll2004r) , g-mode oscillatioris of p rotoneutron stars (PNSs) (lOtt et al.ll2006l ) 
and SASI JKotake et al.lEoOTi : iMarek et al]l2009h . 



Among them, we focused on the GWs o riginated frorn the a sphericities produced by the 



standing accretion shock instability (SASI) (IKotake et al.ll2007l ). Here SASI, becoming very 



shock with 



Scheck et al.l (2004) 



( 120071 ) : llwakami et al 



1 


Blondin et al. 


( 


2003 


): 


Blondin & Mezzacappa 



(l2008al Jbl) and references therein). Based on the two-dimensional (2D) 



simulations, whic h demonstrate th e neutr i no-driven explo s ions a ided by SASI using the light 



bulb scheme (see IJanka &: Miiller! (119961 ): lOhnishi et al.l (120061 ) for details), it was pointed 



out that the GW amplitudes from anisotropic neutrino emission increase almost monoton- 
ically with time, which are dominant over the ones from matter motions, and that such 
signals may be vis ible to next-generation detectors for a Galactic source. More recently, 
Marek et al.l (120091 ) analyzed the GW emission based on their long-term 2D Boltzmann 
(ray-by-ray ] simulations, which s eem very close to produce the SASI-aided neutrino-driven 
explosions (IMarek &: Jankall2009l ). They also found that the GWs from neutri nos with con- 



tinuou sly growing amplitudes (but with the different sign of the amplitudes in IKotake et al. 
( 120071 )). are dominant over the ones from matter motions. They proposed that the third- 



generation class detectors such as the Einstein Telescope are required for detecting the GW 
signals with a good signal-to-noise ratio. 

Except for the acoustic mechanisrn . ( Ott et al. 20061) and the magnetohydrodynamic 
mecha n ism (e.g.,lYamada fc SawaJ (2004 ): 



(120061 ): iBurrows et all (120071 ): iTakiwaki et al 



Kotake et al.l (120041 ): I Ardeli an et all (120051 ) : IShibata et al. 



( 2OO9I ) and references therein), both of which 



-4- 



produce strong mass-quadrupole GWs, all the studies in the postboun ce phase me ntioned 
above, rely basically on the conventional neutrino-heating mechanism (jBethdll99Cll ). They 
agree that the GW amplitudes from anisotropic neutrino emission are dominant over the 
ones from mass motions. This means that accurate estimation of the neutrino anisotropy is 
indispensable for understanding the gravitational radiation from core-collapse supernovae, 
which requires to estimate precisely the directional dependence of the neutrino intensity 
emitted from the central cores . In the previous s imulations with the light-bulb approxi- 
mations (jMiiller fc Jankal 119971 : iKotake et al.l 120071 ) , there was no way but to estimate the 
angle-dependent neutrino luminosities, assuming that neutrinos are emitted purely radially 
in each lateral bin of the computational polar grid. This ray-by-ray treatment cannot cap- 
ture the neutrino emission in t he lateral directions entirely. The FLD (flux -limited-diffusion) 
schemes (either raulti-energv (Walder et al.ll2005l ) or single-energy FLD (IBurrows fc Hayes 



19961 : iFryerl l2004l : iFryer et al.l |200J)) tend to smooth out the local and glo bal neutrino 



anisotropies due to their diff usion characters. Any ray-by-ray transport schemes (IBuras et al. 



20031 : iMarek &: Jankal 120091 ) . albeit coupled to the Boltzmann transport and thus being one 
of the most sophisticated treatment at present, replaces the 2D neutrino transport with the 
ID transport problem along radial rays in every lateral bin, leading to the overestimation 



of the directional dependence of the neutrino anisotropies (see discussions in iMarek et al. 



( 120091 ) ) . More recently, fully 2D multi- angle Boltzmann transport simulations become prac- 
ticable (lOtt et al.ll2008l ). however too computationally expensive currently to perform the 
simulations, satisfying required number o f the momentum-space angles in order to capture 
the neutrino anisotropies accurately (e.g., Ott ( 20091 )). 



These situations motivate us to propose a method to estimate the neutrino anisotropies 
directly linked to the accurate estimation of the neutrino-originated gravitation al radiation. 



Push e d by the strikin g evidences that support the slow rotation of iron cores (IHeger et al. 



20051 : lOtt et al.ll2006l ). we consider an idealized situation that the neutrino radiation field 



from the PNSs are isotropic. This means that no GWs are assumed to be generated inside 
due to the isotropy. Then we focus on the asphericities outside the PNSs, which are produced 
by the growth of SA SI. To this end, we utilize the 2D SASI models as in our previous paper 
(IKotake et al.ll2007l ). in which the neutrino lumino sities from the PNSs are adjusted to trig- 
ger explosions aided by SASI (lOhnishi et al.ll2006l ). Since the regions outside the PNSs are 
basically optically thin to neutrinos, we solve the transport equations by making use of the 
ray-tracing (long-characteristics) method in a post-processing manner. For the sake of sim- 
plicity, we consider here only neutrino absorptions and emission by free nucleons, dominant 
processes outside the PNS, while the neutrino scattering and the velocity-dependent terms in 
the transport equation are neglected. We then study how the obtained gravit ational wave- 
forms could change from the ones obtained in the ray- by-ray analysis (e.g., iKotake et al. 
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( 120071 )) and discuss their implications. Although the presented scheme is designed to be 
valid for the regions where neutrinos are thin to matter, we hope it to be useful for the GW 
prediction in the first generation of 3D core-collapse supernova simulations that do not solve 
the transport equations as part of the numerical evolution. 

The plan of this paper is as follows. In section [2], first the formalism for calculating grav- 
itational waveforms is described. Thereafter, we describe the ray-tracing method to calculate 
the direction-dependent neutrino luminosities. In Section [21 we shall briefly summarize the 
information how to construct the 2D exploding models, such as about the initial models and 
numerical methods. The main results are shown in Section HI We summarize our results and 
discuss their implications in Section [51 



2. Computing the Gravitational Wave Signatures 

2.1. Formulae for Gravitational Waves from Anisotropic Neutrino Radiation 

To compute the gravitational wavefo rms from aniso tropic neutrino radiation, we follow 
the formalism pioneeringly proposed by lEpsteinl (119781 ). In the following, we present the 
formulae in a useful form, which is applicable also for 3D computations. 

The two polarization states (of -|- and x modes) of GWs fro m anisotropic n e utrino 
radiation satisfying the transverse-traceless conditions are given by iMiiller fc Jankal (119971 ) 
as follows. 



2G /"* f 

h+ = — dt' dn'(l + cos9) cos 20 

C'R Jo Jin 



and 

h- = ^iy^' I dn\l + cos9) sin 20 '^^'f^f\ (2) 

where G is the gravitational constant, c is the speed of light, R is the distance of the source 
to the observer , dl^/ dQ represents the direction-dependent neutrino luminosity emitted per 
unit of solid angle into direction of Q. Variables with dash such as Q' represent the quantities 
of the source coordinate system, while non-dashed ones belong to the the observer coordinate 
system (for the geometrical setup, see Figure [T|). For convenience, we assume that y axis 
coincides with y' axis and that the 2;-axis lies on the {x', z') plane. 

Using the following relations between the two coordinates, 

sin 6 cos = sin 9' cos 0' cos ^ — cos 6' sin ^, (3) 



sin 6 sin = sin 9' sin 0', 



(4) 
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Fig. 1. — Source coordinate system {x',y',z') and observer coordinate system {x,y,z). The 
observer resides at the distant point on the 2;-axis. The viewing angle is denoted by C, 
which is the angle between z and z' axis. The z' axis coincides with the symmetry axis 
of the source, presumably the rotational axis. Central red region illustrates the anisotropic 
neutrino radiation from a core-collapse supernova. 



cos 6 = sin 6' cos (p' sin C, + cos 9' cos ^, (5) 

9 and in equations ([H [2]) are required to be expressed in terms of the angles 9', 0' with 
respect to the source coordinate valuables, and the viewing angle of ^. After some algebra, 
one can readily derive the two modes of the GWs (equations ([H [2])) as follows. 



and 



2G /"* f 

h+ = — I dt' + sin6''cos0'sin^ + cos6''cosO X 

C R Jo Jatt 

(sin 6*' COS 0' cos^ - cos6'' sin^^ - sin^ 9' sin^ 0' dl^{Q',t') 
(sin 9' cos 0' cos ^ — cos 9' sin + sin^ 9' sin^ 0' dO,' 



AG /"* f 

= —m / / (^^^'(1 + sin^^'cos0'sin^ + cos6''cos^) X 

C^R Jo Ji^ 

sin^'sin0'(sin^'cos0'cos^ — cos ^' sin ^) dl^{Q',t') 
(sin 9' cos 0' cos ^ — cos 9' sin C,y + sin^ 9' sin^ 0' d^l' 



(6) 



(7) 



which will be useful in computing the GW signals also for 3D computations. It is noted that 
the sum of the squared amplitudes + j/^xP is invariant under the rotation about the 
z-axis. To maximize the GW amplitudes in our 2D axisymmetric case, we as sume that the 



observer is situated along the direction of the equatorial plane (^ = 7r/2) as in lKotake et al. 



- 7- 




Fig. 2. — An example of ray-tracing of neutrinos for estimating dl^{fl) / dQ towards a given 
direction of fl. The central region colored by red represents the surface of the protoneutron 
star (PNS) (located at 50 km in radius), which is the inner boundary of our computation. 
The color-scale on the rays shows the logarithmic values of f^,, the neutrino occupation 
probability (for e^, = 14 MeV here), which is calculated by the line integral along each ray 
(equation (fTTj) ). For this snapshot taken from our 2D simulations, the higher values of f,, 
are seen to come just in front of the PNS (yellow on the rays), while fi, becomes smaller 
in the distant regions from the PNS. Doing ray-tracing calculations along every direction, 
dl^{fl)/dQ can be estimated through equation ([H]). 



( 120071 ). Then the only non- vanishing component is, 



2G f f . ^,^cos''9'-sm'e'sm'<p'dU9',t' 

h,=hy = —— dt (iiZ (1 + smb' cos I ^ 



c^RJo J 4^ ^ ^ cos2 ^' + sin^ ^' sin^ 0' dQ' 



d^RJa Jo ' ' dn' 

Here the subscripts of ('^) and indicates that the observer is situated in the equatorial 
plane and that the GWs are originated from neutrinos. $(^') depends on the angle measured 
from the symmetry axis {6') 

^{e') = TTsme'{-i + 2\cose'\). (9) 
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As given in Figure 1 of iKotake et al.l (120071 ). this function has positive values in the north 
polar cap for < 6' < 60° and in the south polar cap for 120° < 9' < 180°, but becomes 
negative values between 60° < 9' < 120°. In order to perform numerically the angular 
integration in equation (8) with high numerical accuracy, we recommend to perform it based 
on a Gaussian quadrature, because the function of $ of equation ([9]) is not smooth near 
poles and equator. For estimating hy (equation (8)), we are yet to determine the directional 
dependent neutrino luminosities of dly/dVL, which we will estimate by the ray-tracing method 
in the next section. 

Here it is noted that GWs generated by neutrinos are distinct from the ones from matter 
dynamics, because the former has the memory effect, which means that the gravitational 
amplitude jumps from zero to a nonvanishing value and it keeps the non- vanishing value even 
after the energy source of gravitational waves disappeared (see iBraginskii fc Thornd (119871 ) 
for details). In equations (fT|2|) . this nature can be directly seen as the time-integral. Other 
astroph ysical emitters of GWs with memory have be en elaborately s tudied for gamma- ray 
bursts ( Hiramatsu et al. 2005 : Suwa fc Murase 2009), Pop III stars ( Suwa et al. 2007al lbl). 
and inspiralling compact objects ( FavataT2008 ) (see references therein). 



As for the gravitational waves of the quadrupo le radiation of mas s moti ons, we em- 
ploy the stress formula (see, e.g., equation (12) in iMonchmeyer et al.l (Il99ll )). In using 
the formula, we consider the self-gravity of matter in the accretion flow. In the following 
computations, we assume that the source is assumed to be located at our galactic center 
[R = 10 kpc) and also that the observer is situated in the direction of the equatorial plane. 



2.2. Ray- Tracing Calculations of Anisotropic Neutrino Luminosities 

Now we proceed to determine dlu/dVL in equation (8), the directional dependent neutrino 
luminosities through the ray-tracing method. Note in the following equations that we change 
dashed variables to non-dashed ones for the sake of simplicity. 

In the ray-tracing approach, we consider transfer along the ray specified by a constant 
impact parameter p. The coordinate along p is called s, satisfying 

r = (/ + ,2)1/2^ (^0) 

where r is the radial coordinate. To estimate dly{Vt) / dVt along a given direction of fi, we 
draw rays of neutrinos as shown in Figure [2J As will be discussed soon later, in order to get 
the numerical convergence of dljj{Q,) / dVL, we need to set 45,000 rays for each direction, which 
consists of 500 x 90 rays, where the former is for the impact parameters covering from the 
inner- {p^■^^ = 50 km) to the outer- boundary (pout = 2000 km) of the computational domain 
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and the latter is for covering the circumference (e.g., 27r) of the concentric circles on the 
plane perpendicular to the rays (see Figure [2]). In the axisymmetric case here, we perform 
the ray-tracing calculations 60 times to cover the entire sphere, which is the number of the 
mesh points for the polar direction (e.g., section 3). 

The transfer equation of the neutrino occupation probability f,^{eu,p,s) for a given 
neutrino energy along each ray is given by. 



j{e^,p, s)(l - fu{e^,p, s)) - 



fu{eu,P,s) 



(11) 



as A 

where j and A is the emissivity and abs o rptivi t y via neutrino absorp tions and emission by 
free nucleons {ue + n^ e~ + p) (jBruennl f 19851 ): lOhnishi et al.l (l2006l )). which are dominant 
processes outside the PNSs. The optical depth for those reactions are estimated by t,^ = 
1/A. For the sake of simplicity, the neutrino scattering and the velocity-dependent terms 
in the transport equation are neglected here. Although z/g's are taken into account the 
hydro simulations, we focus only on electron-type neutrinos here for simplicity, since they 
dominantl y contribute to the re sulting GWs outside the PNSs as will be discussed in section 
4 (see also iKotake et al.l (120071 )). We use 16 neutrino energy bins which is logarithmically 
uniform and covers 0.9 - 110 MeV. Along each ray, f^, is transferred by the line integral. 
When the line integral starts from the surface on the PNS (the lines coming from the PNS 
in Figure [2]), we set the initial value of 

1 1 

1 + cxp{e^ / k^T^) 47r' 



fie.) 



(12) 



assuming that the neutrino distribution function at the surface is approximated by the Fermi- 
Dirac distribution with a vanishing chemical potential. Here the neutrino temperature is set 
to be constant near T.^ = 4 MeV, whose values change slightly depending on the input 
neutrino luminosity. Note that these values are consta nt in time for each model. This is 
necessary to realize the steady unperturbed states (e.g.. lOhnishi et al.l (120061 )). For the rays 



that do not hit the PNS, we start the line integral from the outer most boundary antipodal 
to the line of sight, where f,^ is essentially zero. 

By a post process, we perform the line integral up to the outer-most boundary for each 
hydro-timestep. The time sampling for the postprocessing is about ~ 1 ms, which is sufficient 
here because the waveforms from neutrinos show much slower temporal variation (> 50 ms) 
as we will discuss in section 4. Since the regions outside the PNSs are basically thin to 
neutrinos (e.g., section 4.1), the transport equation is suited to be solved by making use of 
the ray-tracing (long-characteristics) method. For all the rays, our treatment in equation 
(11) is identical to neglect the time-retardation, which is equal to the light traveling time 
of neutrinos from the PNS to the outer boundary. So the difference of the light traveling 
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time comes from the difference in the position of the neutrino emitting regions, which are 
close to the PNS surface of ~ 50 km. Divided by the speed of hght, the difference of the 
time-retardation, is less than the order of 0.1 ms, which is negligible for our computation, 
because the hydrodynamical timescale is longer than ~ 10 ms. 

With f{eu, p, Sout), which is obtained by the line integral up to the outer-most boundary, 
the neutrino energy fluxes along a specified direction of f2 can be estimated, 

/(e^,p, Sout) ■ (ce^) ■ TTT^T^- (13) 



By summing up the energy fluxes with the weight of the area in the plane perpendicular to 
the rays, we can find dl^/dVL along a specified direction O, 

dS = dp 2-Kp ,^ ■ (14) 



dn I dndS /„. d^ldS 

Repeating the above procedures, dly{Vt)/dVL can be estimated for all the directions. 

To verify the newly developed numerical code, we calculate dlv{Vt)/dVL in a spherical 
medium as a test calculation. The left panel of Figure [3] shows how the deviation of dlyiVt) / dQ 
from its (polar) angle- averaged value, indicated by 6{dl^/dQ) = dl^/dQ — {dl^/dQ)a_^erage in 
the figure, changes with the number of the neutrino rays. With increasing the number of the 
rays, the deviation is shown to be smaller especially near the polar regions (near 0°, 180°), 
while the deviation near the equator (near 90°) is nearly converged, but with relatively 
larger errors than for the other direction. This is because the angular grid size near the 
equatorial belt becomes larger for the spherical coordinate system taken in our simulation. 
Right panel of Figure [3] shows the resulting GW amplitudes for the different number of the 
rays. It is noted that the amplitudes vanish formally for the isotropic configuration. From 
the restriction of the computational time, we choose to set 45,000 rays for our actual ray- 
tracing calculation (green circle in the right panel). With this choice, the numerical errors 
are suppressed to be an order of 10"^'', which is typically 4-5 orders of magnitudes smaller 
than typical GW amplitudes obtained in the models computed here. We have also checked 
the numerical convergence for a deformed matter configuration (see Figure H]) . 



3. Construction of 2D Exploding Models 



As in iKotake et al.l (120071 ) , we assume that the explosion is powered by the neutrino- 
energy deposition between the PNS and the standi ng accretion sho c k, pro bably one of the 
most promising way to blow up massive stars (e.g., iMarek &: Jankal (120091 )). To trigger the 
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Fig. 3. — Numerical tests, in which the ray-tracing calculation is done for a spherical medium 
using data at t = 294 ms for model A. Left panel shows the deviation of dl^{Q)/dQ from 
its angle-averaged value as a function of the angle measured from the polar axis for different 
number of the neutrino rays (:A^ray) (see text for more details). Right panel shows the 
resulting GW amplitudes, which should formally vanish. The green circle represents for 
Nr»^, = 45, 000, which we employ for our actual ray-tracing calculation. Note in the right 



ray 



panel that the supernova is assumed to be located at a distance of 10 kpc. 



explosions, we empl oy the so-called light-bulb approximation (see IJanka &: Miilleii (Il996l ): 
Ohnishi et al.l (120061 ) for details) and adjust the neutrino luminosities from the PNSs, which 
was found t o work well iii 2D, capturi ng the main features obtained by more realistic simu- 



lations (e.g.. IScheck et al.l (l2004l . l2006l )). It should be noted again that we focus on the GWs 



generated by the anisotropies associated with the growth of SASI, whic h develops outside 
the PNS. Therefore, the oscill ations of the PN Ss (IBurrows et al.l l2006l ) and the resulting 
efficient gravitational emission Jott et ahlbooeh or the enhanced neutrino emission near the 



equator of the PNS observed in iMarek et al.l ( 120091 ) . cannot be treated in principle here. 

I992I) and 



Our hydrodynamic code is based on the ZEUS (iStone &: Norman 



maior 



modifications for th e supernova simulations are already described in iKotake et al.l (120031 ): 
Ohnishi et al.l (120061 ). The computational grid is comprised of 300 logarithmically spaced ra- 
dial zones to cover from the absorbing inner boundary of ~ 50 km to the oute r boundary of 
2000 km, and 60 polar (6) uniform mesh points (see for the resolution tests in iKotake et al.l 
J2007h l 



The initial condi tion is a spherically symme tric steady accretion fiow through a stalled 
sh ock wave to a PNS ([Yamasaki fc Yamada 2005), which is also provided in the same manner 



of lOhnishi et al. 

pin = 10" g cm 



J2006h . In constructing the initial conditions, we assume a fixed density 
at the inner boundary. And the initial mass accretion rates and the initial 
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1e-21 



100 1000 10000 

Number of Neutrino Rays 



100000 



Fig. 4. — Same as the right panel of Figure [3] but for a deformed medium. Here we artificially 
impose a global quadrupole deformation {i = 2) to the matter distribution. It shows the 
numerical convergence of the GW amplitudes for A^^ay = 45, 000. 



mass of the central object are set to be M = 1 Mq s ^ and M^^ = 1.4 Mq, respectively. 



As in iKotake et al.l (120071 ) . the self-gravity of matter in the accretion flow is ignored. To 
induce non-spherical instability, we add random velocity perturbations of less than 1 % of 
the unperturbed velocity. At the outer boundary, we adopt the fixed boundary condition 
consistent with the initial condition. On the other hand, the absorbing boundary is used at 
the inner boundary. We assume that the neutrino flux from the PNSs can be approximated 
by black-body emission. 



Pushed by the evidences that the SASI-induced explosions are f avorable for exp 



the observed qua ntities such as the synthesized elements of SN1987A (IKifonidis et al 



the pu lsar kicks (jScheck et al.ll2004j ) and spins (IBlondin fc Mezzacappa 



2007 



ainm g 



2006), 



Iwakami et al. 



2008bl ). we mainly focus on the gravitational radiation in the models tuned to produce 
explosions. By changing the electron neutrino luminosity at the surface of the PNS in the 
range of L,y^ = 6.4 — 6.8 x 10^^(10 foe) erg s~^, we construct three exploding and one non- 
exploding models (see Tabled]). Except for model D, we can observe the continuous increase 
of the average shock radius with the growth of SASI, reaching the outer boundary of the 
computational domain with the explosion energy of ~ 10^^ erg. Until this moment, we run 
the simulations (e.g.. At in the table), while in model D, we terminated the simulation at 
about 800 ms, not seeing the increase of the shock radius. 
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4. Results 



The left panel of Figure [5] shows the GW amplitudes contributed from anisotropic neu- 
trino emission for different luminosity models. Comparing the right panel, which shows 
the total amplitudes (neutrino + matter), it can be seen that the gross structures of the 
waveforms are predominantly determined by the neutrino-originated GWs with the slower 
temporal variations (> 50 ms), to which the GWs from matter motions with rapid temporal 
variations (< 10 ms) are superimposed. 

On the other hand, the waveforms shown in the left panel of Figure [6] are estimated 
simply by sur aming up the local neutrino cooling rates outside the PNSs with the ray-by-ray 
assumptions (IKotake et al.ll2007l ) (see also Figure 2 in the paper). They are so much different 
from the ones obtained here. In contrast to the monotonic increase of the amplitudes with 
time, the waveforms here exhibit more variety, showing large negative growth at some epochs. 
And the GW amplitudes from neutrinos become more than one-order-of magnitude smaller 
than the previous estimation (compare Figures O and the left panel of E] noting the different 
vertical scales). On the other hand, similarities between them are that the waveforms from 
neutrinos have generally a positively growing feature with time, and also that the electron- 
type neutrinos dominantly contribute to the wave a mplitudes than for th e anti-electron type 
neutrinos (right panel of Figure El see also Figure 3 in lKotake et al.l (120071 )). In the following, 
we first analyze this positive trend based on the ray-tracing calculations. And in section 14. 2[ 
we clarify the reason for the negative growth, contributing to make the amplitudes much 
smaller than previously estimated. 



Table 1. Model Summary 



Model L^, (1052 erg/s) At (ms) h^,fi„ (IQ-^^) |/itot,max| (IQ-^^) Eaw,^ (IQ-i^Mqc^) 



A 6.8 509 8.7 7.7 0.44 

B 6.7 570 2.2 9.1 1.32 

C 6.6 740 6.1 8.0 1.39 

D 6.4 800 4.8 6.1 0.49 



Note. — L^^ denotes the input luminosity. At represents the simulation time. /i^,fln and /itot,max 
represents the amplitudes of the neutrino-originated GWs at the end of the simulations and the maximum 
amplitudes (neutrino + matter) during the simulation time. i?GW,i^ is the radiated energy in the form of 
the neutrino GWs in unit of Mqc^. Note that the supernova is assumed to be located at a distance of 10 
kpc. 
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Fig. 5. — Gravitational waveforms only from anisotropic neutrino emission (left) and from 
the sum of neutrinos and matter motions (right). The time is measured from the epoch 
when the neutrino luminosity is injected from the surface of the neutrino sphere. In all the 
computed models, SASI enters to the non-linear regime at about 100 ms, simultaneously 
making the amplitudes deviate from zero. Note that the supernova is assumed to be located 
at the distance of 10 kpc. 



4.1. Positively growing features 

To see the positively growing features clearly, we choose model A as a reference. The 
left panel of Figure [7] shows dlu/dVL in the vicinity of the north pole {6 = 0), the equator 
{9 = vr/2), and the south pole {9 = it). In the following, we focus on the two epochs of t 
= 294 and 370 ms, when the neutrino emission in the south pole dominate over the ones 
in the north pole and vice versa. The right panel of Figure [3 shows that the dominance of 
the neutrino emission in the north and south poles are closely anti-correlated. This is the 
consequence of the low- mode nature of SASI, here of £ = 1. In fact, the right panel of Figure 
E] shows that at 294 ms, the blob encompassing the regions inside the stalled shock is moving 
from the southern to the northern hemisphere, leading to the compression of the matter in 
the south hemisphere, which is vice versa at 370 ms (right panel). 

Figure M shows various properties obtained by the ray-tracing calculation for t = 294 
ms. Top, middle, and bottom panel, is the case seen from the northern hemisphere, the 
equator, and the south hemisphere, respectively. Comparing the right panels, we can see 
that the rays with highest fi^, come from the high-temperature regions near the vicinity of 
the south pole (bottom right panel). As mentioned above, this is due to the compression of 
matter near the south pole, produced by the sloshing SASI, moving from the south to the 
northern hemisphere. 
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Time [ms] Time [ms] 



Fig. 6. — Gravita tional waveforms fro m anisotropic neutrino emission estimated by the ray- 
by-ray analysis of iKotake et al.l (120071 ) . showing a monotonia increase of the amplitudes with 
time (left). Right panel is the same as the left panel of Figure |5] for model A but with the 
contribution of anti-electron type neutrinos [ug) to the waveform (indicated by "Total"), 
showing that the dominant contribution to the GW amplitudes comes from electron-type 
neutrinos (z/g). The supernova is assumed to be located at the distance of 10 kpc. 



In the left panels of Figure [9l it is noted that due to the axisymmetry, the distributions 
of the neutrino energy fluxes seen from the north (top left) and south hemispheres (bottom 
left) have the circumferential structures. The images of the energy fluxes may seem like an 
annular eclipse, which we explain as follows. Here for clarity, let's consider two rays of A 
and B. The ray A is coming a point from {X = 0,Y = 0) in the middle left panel, which 
in the middle right panel, is equal to the point along the line perpendicularly threading the 
center of the colored plane, but is positioned at 50 km away from the plane to the observer 
(:right direction in the panel). The ray B is coming from a point from (X = 0,Y ^ —50km) 
in the middle left panel, which is equal to the point, near the southern point on the colored 
plane in the middle right. The ray B travels to the observer, experiencing longer the high- 
temperature regions near the surface of the PNS. This makes the eclipse-like shining near 
the edge of the PNS. The left panel of Figure [TU] is an enlargement of the bottom right panel 
of Figure [9] near the PNS, which clearly shows the rim-shining near X = 50 km (compare 
the right-hand side with the bottom left panel of Figure [9]). The right-hand side in the right 
panel of Figure [10] depicts the optical depth of neutrinos at e^, = 14 MeV. It is seen that the 
high values of fi, come from the higher temperature regions in the vicinity of the southern 
part of the PNS (left panel), where are thin to the neutrinos with the optical depth being 
smaller than 2/3 (right panel). Such an optically thinness is favorable for the ray-tracing 
calculation here as discussed in section 12.21 
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Fig. 7. — The directional dependent neutrino luminosity: dly/dTt for model A, in the vicinity 
of the north pole, the equator, and the south pole (left panel), and their differences from the 
equator (right panel). Vertical lines represent the epochs of t = 294 and 370 ms, respectively 
(see text for details). 



Among the left panels of Figure [HI the highest values of the local neutrino energy fluxes 
{:dlu/{dQdS)) come from the brightly shining southern poles seen from the equator (middle 
panel). However by summing up them with the weight of the area (e.g., dS in equation 
f|T^ ). dl^/dQ are shown to become largest seen from the southern hemisphere (the red line 
in the left panel of Figure [TT]). The similar argument is true for t = 370 ms (Figure [T2]) . 
while the neutrino luminosities seen from the northern hemisphere are higher than the ones 
from the southern hemisphere (the green line in the left panel of Figure [TT]) . Noting again 
that $(6' ) in equation ([8]) are positive near the north and south polar caps, the dominance 
of the anisotropic neutrino luminosities in the vicinity near the north and south poles, makes 
the positively growing feature in the resulting GWs. 



Figure [13] shows time evolution of a neutrino anisotropy parameter defined in lMiiller &: Janka 



(119971), 

where Lj,(t) is the total neutrino luminosity L,y{t) = j^^dVt'dl^{Vt' ,t')/dVL' . This quantity is 
useful to see how the temporal changes of the neutrino anisotropy have impacts on the GW 
amplitudes. Either for models A or B, it can be seen that a keeps positive value with time 
in the later phase (> 400 ms). This coincides with the epoch when the low-modes explosion 
is triggered by SASI along the symmetry axis, which is also helpful for understanding the 
reason of the positive growth (compare the left panel of Figure [5]). 
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Fig. 8. — Temperature- (the left-half of each panel) and density- (the right-half) distributions 
in the meridian section for model A. The insert of each panel shows the neutrino GW 
amplitudes, in which the green point indicates the time of the snapshot. The central region 
colored by white (50 km in radius) represents the inner boundary of our computation. 



4.2. Negatively growing features 



Now we move on to discuss large negative growth of the GW amplitudes seen in the 
waveforms (e.g., models B and C in Figure El), wh ich is newly found by the ray-tracing 



calculation (compare Figure 2 of lKotake et al.l (120071 )). Here we choose a snapshot at t = 330 



ms of model B as a reference, which shows the steep negative growth (see the insert of Figure 



The middle panel of Figure [T^ shows that the neutrino emission seen from the equatorial 
plane are much stronger than the ones seen from the northern/southern hemispheres. The 
bright belt-like structure in the middle panel is from the high-temperature regions in the 
northern hemisphere, which are by chance surrounded by relatively low density (see cusp-like 
structure in the northern hemisphere of Figure [T^ and [TBI) . In this case, the area- weighted 
energy fluxes become larger seen from the equatorial plane (see right panel of Figure [TT]) . In 
combination with the negative values of ^{O') in equation ([8]) in the vicinity of the equatorial 
belts, this makes the negative growth in the GW amplitudes. Large negative amplitudes seen 
for some other epochs in other model such as model C (left panel of Figure [5]) are also from 
the same reason. Such a feature is genuine outcome of the neutrino emission in the lateral 
direction, which is able to be captured correctly by the ray-tracing calculation. 
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It is noted that the appearance of the negative growth has no systematic dependence 
of the input luminosities. In fact, as seen from Figure [5], the negative growth is observed for 
the intermediate luminosities models (models B and C), but not for the highest (model A) 
and smallest luminosity models (model D) (see also |/itot,max| in Table 1). This should reflect 
the nature of the SASI which grows chaotically and non-locally. Albeit with the negative 
growth, our results suggest that the positively growing features dominate over the negatively 
ones for the 2D models (see hi,^^ in Table [1]). This is due to the axial symmetry, along which 
the SASI develops preferentially and the resulting anisotropics become larger. 

As mentioned earlier, the neutrino GWs become more than one-order- of magnitude 
smalle r than the previous estimation (compare -Egw, u in Tableland the one in lKotake et al. 



( 120071 )). This stems not only from the incursion of the negative contributions but also from 
the appropriate estimation of the neutrino absorptions made possible by the ray-tracing 
method. Previously the neutrino lumino sity was estimated simply by summing up the local 



neutrino cooling rates outside the PNSs (IKotake et al.ll2007l ). which fails to take into account 



the neutrino absorption correctly (A in equation (fTTl) ). These two factors make the ampli- 
tudes much smaller than the previous estimation. As a result, the neutrino GWs, albeit 
dominant over the matter GWs in the lower frequencies below ~ 10 Hz (Figure [T7|) . become 



very difficult to be detected for ground-basec 
the seismic noises at such lower frequencies 


detectors whose sensitivity is limited mainly by 


Ando & the TAMA collaboration 


2002; 


Thorne 


1995; 


Weinstein 


2002; 


LCGT Collaboration 


1999). 



On the other hand, the GWs from matter motions seem marginally within the detection 
limits of the currently running detector of the first LIGO and the detection seems more 
feasible for the detectors in the next generation such as LCGT and the advanced LIGO for a 
Galactic supernova. The spectra of the matter GWs have double peaks namely near 100 Hz 
and 1 kHz. While the latter comes from the rapidly varying local hydrodynamical instabilities 
with milliseconds timescales, the former is associated with the longer-ter r n ove rturns of 0(10) 
ms induced by ^ = 2 mode of SASI (see e.g.. Figure 5 in iKotake et al.l (120071 )). These gross 
properties in the GW spectra are common to the other luminosity models. Thus the peak in 
the spectra near ~ lOOHz is found to be a characteristic feature obtained in the 2D models 
computed here. 



5. Summary and Discussion 

We proposed a ray-tracing method to estimate gravitational waves (GWs) generated by 
anisotropic neutrino emission in supernova cores. To calculate the waveforms, we derived the 
GW formulae in a useful form, which are applicable also for 3D computations. Pushed by the 
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studies supporting the slow rotation prior to core-collapse, we considered an idealized situa- 
tion that the neutrino radiation field from the protoneutron stars (PNSs) are isotropic. Then 
we focused on the asphericities outside the protoneutron stars, which are produced by the 
growth of standing accretion shock instability (SASI). Since the regions outside the PNSs are 
basically thin to neutrinos, we solve the transport equations by making use of the ray-tracing 
method in a post-processing manner. For simplicity, neutrino absorption and emission by 
free nucleons, dominant processes outside the PNSs, were only taken into account, while 
the neutrino scattering and the velocity-dependent terms in the transport equations were 
neglected. Based on the two-dimensional models, which mimic SASI-aided neutrino heating 
explosions, we estimated the neutrino anisotropies by the ray-tracing method and calculated 
the resulting GWs. 

Our results show that the waveforms from neutrinos exhibit more variety in contrast to 



the ones previously estimated by ray-by-ray analysis (e.g.. iKotake et al.l (120071 )). In addition 
to the positively growing feature, which was predicted to determine the total wave amplitudes 
predominantly, the waveforms show large negative growth for some epochs during the growth 
of SASI. Such a feature is a genuine outcome of the neutrino emission in lateral directions, 
which can be captured correctly by the ray-tracing calculation. Reflecting the nature of 
SASI which grows chaotically with time, little systematic dependence of the input neutrino 
luminosities on the maximum amplitudes and on the radiated GW energies are found. Due to 
the negative contributions and the neutrino absorptions appropriately taken into account by 
the ray-tracing method, the GWs from neutrinos become more than one-order-of magnitude 
smaller than the previous estimation, making their detections very hard for a galactic source. 
On the other hand, we point out that the gravitational-wave spectrum from matter motions 
have its peak near ~ 100 Hz, reflecting the growth of £ = 2 mode of SASI with timescales of 
0(10) ms. Such a feature is found to be generic among 2D models computed here. 

It should be noted that the approximations taken in the simulation, such as the excision 
inside the PNS with its flxed inner boundary and the light bulb approach with the isentropic 
luminosity constant with time, are the very flrst step to model the dynamics of the neutrino- 
heating explosion aided by SASI and study the resulting GWs. As already mentioned, the 
excision of the central regions inside PNSs ma y hinder the efficient gravitational emission 



of the oscillating neutron star (jOtt et al.ll2006l ). It is recently reported t hat the waveforms 



obtained in the 2D ray-by-ray Boltzmann simulations (IMarek et al.ll2009l ) show a negatively 
growing feature. At first glance, our results may seem to contradict with theirs. They 
explained that the negative growth is due to the enhanced neutrino emission of muon and 
tau neutrinos near the equator of the PNS observed in their simulations. As we mentioned 
repeatedly, such feature inside the PNSs, cannot be captured here in principle. Outside the 
PNSs, however, the stronger neutrino emission of Vf, and Vf, along the polar regions are also 
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seen in their simulations (see their bottom right panel of Figure 11 outside of 30 ~ 40 km, 
which is the surface of the PNS guessed from their Figure 1), which we can say, consistent 
with our results. As illuminated by this study, the elaborate estimation of the neutrino 
anisotropy is crucial to understand the gravitational radiation in supernova cores. This 
naturally needs a full transport simulation coupled to 3D hydrodynamics encompassing the 
whole iron core, which is beyond the scope of this study and a grand challenge for all the 
supernova modellers. 

In addition, one more major deficit is the axial symmetry assumed in the present 2D 
simulations. In 3D, the pronounced dominance of / = 1, 2 along the symmetry axis, which 
is a coordinate singularity in the 2D computations, may become weaker, owing to the addi- 
tional spatial degree of freedom in the azimuthal direction. In the 3D case, we expect that 
the amplitudes becorn e smalle r owing to the reduced anisotropy along the symmetry axis 
(llwakami et aP boOSal lbl. boooh . Thus the amplitudes calculated in this study could be an 
upper bound, in which the maximal anisotropy generated by SASI could be achieved. With 
the ray-tracing method and the GW formulae derived here, we ar e now able to investigate 



the properties of the cross mode of GWs (e.g.. iKotake et al.l (120091)). We think it important 



because they are of genuine 3D origin, which in combination of the plus modes, might give 
us son ae hints about the e xplosion asymmetry. Magnetic eff ects not only on the growth of 
SASI (lEndeve et al.ll2008l ). but also on the GW waveforms (IKotake et al.l 120041 ) should be 
interesting. This study is an appetizer before the forthcoming 3D studies to clarify those 
aspects, which will be presented elsewhere soon. 
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Fig. 9. — Various properties obtained by the ray-tracing calculation at t = 294 ms for model 
A. Right panels show the rays (colored by the value of f,^ for e^, = 14 MeV), which are 
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Fig. 10. — Left panel is an enlargement of the bottom right panel of Figure [9] near the 
central PNS (white circle). Right panel depicts the optical depth of neutrinos at = 14 
MeV (right-hand side) with the density distributions (left-hand side). The opacity sources 
are taken to be the neutrino emission/absorption as mentioned in section [221 It can be seen 
that high values of fi, come from the high temperature regions in the vicinity of the PNS, 
where are thin to neutrinos with the optical depth being smaller than 2/3. 
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Fig. 11. — The directional dependent neutrino luminosities of dl^/dQ as a function of the 
polar angle for model A (left panel) at t = 294 and t = 370 ms, and model B (right panel) 
at t = 330 ms. 
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Fig. 12. — Neutrino energy fluxes of dl^ / {dVLdS) (equation (fT3l) ) of model A at t = 370 ms, 
seen from the northern hemisphere (left), the equator (middle), and the southern hemisphere 
(right), respectively. 
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Fig. 13. — Time evolution of the neutrino anisotropy parameter (:a in equation ( ITSi) ) for 
models A and B. a keeps positive value with time in the later phase (> 400 ms) when the 
low-modes explosion is triggered by SASI along the symmetry axis. 
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Fig. 14. — Same as Figure [8] but for model B at t = 330 ms. 
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Fig. 15. — Same as Figure fT2] but for model B at t =330 ms, showing the stronger emission 
seen from the equator (middle), as opposed to the cases of Figures |9] and [T2l 
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Fig. 16.— 



Same as Figure [9] (middle right) but for model B at t =330 ms. 
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Fig. 17. — Spectral distributions of GWs from matter motions (Matter) and anisotropic 
neutrino emission (Neutrino) for r aodel A with the ex pected detectio n hmits of TAMA 



f Ando fc the TAMA collaboration! 120021 ) . first LIGO (IThornel Il995l ). advanced LIGO 



(LCGT Collaboration 
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in 
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^char is the characteristic gravitational wave strain defined 



distance of 10 kpc. 
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